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Abstract. An overview of current status and prospects of the development of 
quantum computer hardware based on inorganic crystals doped with rare earth ions 
is presented. Major parts of the experimental work in this area has been done in two 
places, Canberra, Australia and Lund, Sweden, and the present description follows 
more closely the Lund work. Techniques will be described that include optimal 
filtering of the initially inhomogeneously broadened profile down to well separated 
and narrow ensembles, as well as the use of advanced pulse-shaping in order to achieve 
robust arbitrary singlc-qubit operations with fidelities above 90%, as characterized by 
quantum state tomography. It is expected that full scalability of these systems will 
require the ability to determine the state of single rare-earth ions. It has been proposed 
that this can be done using special readout ions doped into the crystal and an update 
is given on the work to find and characterize such ions. Finally, a few aspects on 
possibilities for remote entanglement of ions in separate rare-earth-ion-doped crystals 
arc considered. 



PACS numbers: 3.65.Wj, 03.67.Hk, 42.50.Md 



1. Introduction 



Many solid state quantum computing systems rely on a top down approach with carefully 
crafted systems, where parameters are designed during the manufacturing or growth 
process, for instance, superconducting qubits [H [2], quantum dots j3] and the Kane 
model jl]. At the other end of the scale there are bottom up approaches, where 
useful hardware can be extracted out of random materials which intrinsically have very 
favorable properties, such as nitrogen vacancy (NV)-centers [5] and rare-earth-ion-doped 
inorganic crystals [6]. In this paper we will discuss development of quantum computer 
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hardware based on inorganic crystals doped with rare earth ions. Such crystals are 
used extensively in current technology, e.g. as active laser materials or as scintillating 
materials in X-ray detectors. They are manufactured in single crystal growth processes 
where a dopant material is added during the growth such that the dopant ions replaces 
host ions in the regular crystal lattice at concentrations that may range between 0.01% 
and 10% depending on application. Although the average concentration is defined 
by the dopant concentration in the crystal melt, the relative positions of the dopant 
ions is largely random. However, randomness and less precise control of construction 
and growth parameters can be compensated for by using suitable quantum computing 
schemes and by developing robust pulses enabling high fidelity qubit operations. 

The next section (Section [2]) will describe basic properties of rare-earth- ion-doped 
crystals relevant for quantum information. This includes coherence times which can 
be several seconds [7] and ion-ion interactions as a mechanism for gate operations. 
For a laser beam passing through a rare-earth-ion-doped crystal there may be 10 15 
ions within the laser focal region. However, qubits can still be created and prepared 
in well defined states and arbitrary qubit operations, characterized by full quantum 
state tomography, give fidelities, F, in the range 0.91 < F < 0.96 [8]. A brief 
description on the experimental setup, including the pulse-shaping system, will be given 
in Section [31 while Section [4] describes how to efficiently initialize the system into well 
defined qubits, which are inhomogeneous subgroups out of these 10 15 ions. In Section [5] 
we show experimental data on qubits prepared in arbitrary superposition states and 
characterized by quantum state tomography. The high fidelity results are all based on 
qubit operations that are robust against individual differences between the ions [9j [10] . 
We will also present data from the development and use of optimal pulses sequences [TT] 
for operation fidelities > 0.96. Section [6] starts with a description of how to, in contrast 
to the experiments mentioned above, create qubits in these crystals that each consist 
of just a single ion, which offers a better scalability. The single ion qubit scheme [12] 
is based on co-doping the crystals with "readout ions". Such readout ions should have 
a short upper state lifetime in order to enable repeated cycling. One such candidate is 
the 4f to 5d transition in Ce 3+ and experimental data to characterize this transition is 
also presented in Section [6j The last Section before the summary contains a discussion 
on entanglement of remote qubits and further scalability. 

2. Rare-earth-ion-doped crystals 

A major strength of using rare-earth ions as qubits are the long coherence times. In these 
ions, transitions within unfilled electron shells are partially shielded from environmental 
dephasing mechanisms by energetically lower, but outer-lying, filled shells. The present 
results were obtained using the 3 H^ — > l D 2 transition in Pr 3+ :Y 2 Si05 at 606 nm [13]. For 
temperatures below 4 K, the Pr 3+ ions in this material has an electronic state lifetime of 
164 [is, and a coherence time, T 2 , of about 100 |os (somewhat depending on the density 
of excited states in the material). By applying magnetic fields however, this time can 
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be increased, and there are also other rare-earth ions which have even longer coherence 
times, for example Er, where a T 2 of 6.4 ms has been measured [T4] . 

The qubit states are not represented by electronically . „ M tj 



inhomogeneous broadening of the absorption line. For 
Pr 3+ , the inhomogeneous linewidth is about 5 GHz, which can be compared to the 
homogeneous linewidth of about 3 kHz, given by the optical T 2 , mentioned above. 

In the rare-earth quantum computing scheme, the inhomogeneous structure is of 
significant importance, as multiple qubits are addressed in the frequency domain, each 
qubit having a different resonance frequency. From above, the ratio of inhomogeneous 
over homogeneous linewidth is about 10 6 for Pr 3+ :Y 2 SiO"5 , but since ions do not have 
only a single transition frequency, the number of frequency channels are in practice 
much smaller. Fig. [1] shows that the lowest and highest transition frequency for a given 
ion differ by 4.8 + 4.6 + 10.2 + 17.3 = 36.9 MHz. Thus, about 100 individual frequency 
channels are a more realistic value. The Pr 3+ ion sits in a non-centrosymmetric site, and 
thus has a permanent electric dipole moment. The value of this dipole moment depends 
on the electron configuration and the value of the dipole moment of an electronically 
excited state, is therefore different from that of the ground state. Qubit-qubit gates can 
then be executed through dipole-dipole interactions. Since the dipole moment changes 
when the ion becomes excited, exciting an ion will shift the transition frequency of all 
near-lying ions. An ion belonging to one qubit can then be controlled by the excitation 
of a another ion, belonging to another qubit, if they are spatially sufficiently close to 
each other. This effect is used as the control mechanism for conditional gates, and the 
basis for this was demonstrated in |10j . 

3. Experimental setup 



excited states though, but rather by ground state hyperfine 
levels, shown in the level diagram of Pr 3+ in Figure [TJ For 
Pr 3+ :Y 2 SiO"5 , the hyperfine lifetime is about 90 s and the 
coherence time about 500 |os. The hyperfine T 2 can also 
be increased by magnetic fields, and up to 860 ms has 
been demonstrated [15]. In addition, dynamic decoupling 
sequences (bang-bang pulses) can also be applied to further 
increase the coherence time to more than 30 s [7]. 




The rare-earth ions are randomly doped into the host 
crystals, replacing one of the host ions in the lattice. 
This replacement causes changes to the crystal electric 
field which then causes each rare-earth ion to pick up a 
shift in its resonance frequency. This amounts to a large 



Figure 1: Relevant en- 
ergy level structure for 
the Pr 3+ ion transition 
utilized in the present 
work. 



The experimental system is based around a Coherent 699-21 ring dye laser, stabilized to 
1 kHz linewidth, by locking to a spectral hole in a Pr 3+ :Y 2 SiOs crystal [16] . The rare- 
earth crystal used in the experiments is a 0.5 mm thick Pr 3+ :Y 2 SiOs crystal with a Pr 3+ 
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doping concentration of 0.05%, which gives a maximum optical density of about 2 — 3. 
In order to get rid of phonon interactions, the crystal is submerged in liquid helium and 
cooled to 2 K. The available laser power is in the order of 100 mW before the cryostat, 
which translate into maximum Rabi frequencies of about 1-3 MHz, depending on which 
transition is targeted, given a laser focus of 100 [im. 

High fidelity gates requires a high-performance pulse-shaping system, as will be 
further discussed in Section |5j Two acousto-optic modulators (AOMs) are used to 
create the desired pulses. The incoming continuous wave (CW) photons are deflected 
by an RF acoustic wave in the AOMs. By varying the power of the acoustic wave, 
the amplitude of the deflected light is controlled. By virtue of energy and momentum 
conservation, the deflected photons must increase their frequency by an amount equal 
to that of the phonons. Thus, by varying the frequency of the acoustic wave we directly 
control the frequency of the light pulse, and furthermore, by changing the phase of the 
acoustic wave we also gain control over the optical phase of the light, which is necessary 
for controlling the coherent interaction with the ions. For these modulators, the power 
conversion from RF wave to light amplitude is quite non-linear. In order to achieve high 
accuracy of the pulse amplitude envelopes over as much as six orders of magnitude in 
dynamic range, an elaborate calibration system, that also includes the power dependence 
on RF frequency, was setup. The first AOM has a center frequency of 200 MHz and is 
aligned in a double pass configuration, in order to cancel out spatial movement of the 
deflected beam. After both passes the AOM has a total deflection range of 200 MHz. 
The second AOM has a 360 MHz center frequency, and is used only to split a pulse at 
a single frequency into two frequencies, such that an excitation pulse can be resonant 
with the |0) — > |e) and |1) — > |e) transitions at the same time, as further discussed in 
Section 15.11 

4. Qubit initialization 

Starting from the inhomogeneous distribution, we employ advanced optical hole-burning 
sequences in order to extract a useful qubit system. Initially, the ions occupy all three 
ground states, and the absorption at any given frequency within the absorption profile, 
contains contributions from all nine possible transitions (see Figured]). The goal of the 
hole-burning is to initialize a narrow controlled ensemble of ions into the |0)-state, well 
separated in frequency from any other ion. 

At first, a wide spectral hole, henceforth called a pit, is created. The hyperfine 
level splittings set the limit for the maximum frequency range of such a pit. Since the 
ions must be in one of the ground states, it can not be larger than 27.5 MHz, as can 
be seen in Figure [U This maximum width is then further reduced by the excited state 
splitting down to 18.1 MHz. Scanned laser pulses target specific transitions inside this 
18 MHz region, and are repeated in order to create a spectral pit where all ions have 
been removed within an interval of maximum 18 MHz. Such a simple pit is shown in 
Figure |2|l. 
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Figure 2: (color online) a) shows a simple pit created by repeated scanning across an 
interval less than 18 MHz. Panel b) displays an improved pit, where scans on multiple 
intervals have been iteratively repeated in order to shuffle all ions outside the pit as far 
out as possible, thereby creating an optimal pit. Also note that the scales in the two 
experimental figures are different. 



Eventually, qubit operation pulses will be used inside the pit structure, and it then 
becomes important to minimize the interactions between the tails of such pulses and 
the edges of the pit. The pit shown in Figure [2^ displays non-optimal edges, that are 
a bit rounded and also go up to the maximum aL immediately. In order to improve 
on this, a more elaborate pit-burning scheme was introduced. In this scheme, pulses 
that target specific levels of Pr 3+ outside the pit are used in order to burn them even 
further out, as illustrated by the top part of Figure |2b. For an ion having a particular 
frequency relative to the pit, as the ion illustrated in this figure, level |0) and |1) are 
targeted in an optical pumping process that eventually transfers them to their \aux) 
state. After this, ions of other frequency ranges might fall back down inside the pit, 
and the initial pit-burning pulses will have to be repeated iteratively together with the 
secondary pulses to keep the center of the pit empty. A similar procedure is then also 
repeated on the left side of the pit in order to achieve an optimal result, shown in the 
lower part of panel b. A more extensive description of the burning procedure can be 
found in [IT] and a detailed list of the exact pulses used to burn such an optimal pit 
can be found in [IB] . 

After an empty spectral pit has been created, a qubit ensemble can be initialized 
inside the pit by a "burnback" pulse sequence. These pulses can with benefits be coherent 
pulses that transfer ions from outside the pit, say at higher frequencies, first to the 
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Figure 3: (color online) a) shows the rare-earth qubit system initialized to the |0) state, 
visible through the three absorption peaks from this state. Panel b) then shows the 
system after a coherent transfer over to the |1) state. 



excited state, and from there down to one of the ground states in a controlled manner. 
Figure [3^ shows the system initialized to one of the ground states, in this case the |0) 
state. 

5. Qubit gate operations 

The ensemble qubit is now initialized to the |0) state. In order to be able to 
perform high fidelity gate operations, we must now find pulses that compensate for 
the inhomogeneous broadening of the ensemble, which would otherwise cause dephasing 
during the excitation. One type of pulses that are robust against such dephasing are the 
complex hyperbolic secant pulses, or sechyp for short. The amplitude envelope of these 
follow a hyperbolic secant shape, while at the same time the pulse is chirped according to 
a hyperbolic tangent function, centered around the mean ensemble transition frequency. 
Such pulses can efficiently transfer an inhomogeneous ensemble between the two poles 
on a Bloch sphere, as illustrated in Figure Hh. The different lines correspond to ions 
with different atomic resonance frequencies, which can all be seen to converge at the top. 
The robustness of these pulses has been verified experimentally by transfer efficiencies 
of 97% for ensembles inhomogeneously broadened more than 100 linewidths. Figure [3b 
shows the system after such a sechyp transfer from |0) — > |e) immediately followed by 
a transfer from |e) — > |1). In addition to this, the sechyp pulses also have the good 
property that they are robust against fluctuations in laser power. While normal square- 
or gaussian-shaped pulses need to have very well tuned values of laser intensity and 
duration in order to get the correct pulse area, the sechyp pulses simply need to be 
above a certain threshold Rabi frequency, above which they always work with high 
efficiency. This independence of the light intensity can be intuitively perceived from 
the Bloch sphere in Figure Hh, where we see that as the different components approach 
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Figure 4: (color online) a) shows the paths taken by three ions with different resonance 
frequencies, while being transferred to the excited state by a coherent sechyp pulse. 
Panel b) illustrates the dark and the bright states on the Bloch sphere, relevant for 
the gate operations, and panel c) shows a full qubit gate operation by means of two 
consecutive dark state pulses (two-color pulses). Note that the transfer in c) is going 
from and to the bright state. 



their final state, they circle around the top pole, and a higher Rabi frequency would just 
make them circle around the z-axis faster, not take them away from the desired state. 
These pulses were investigated more carefully in [9]. 

5.1. Dark state gates 

The sechyp pulse-shapes solve the problem of inhomogeneous dephasing as well as offer 
a robustness with respect to fluctuating experimental parameters, but can only perform 
pole-to-pole transfers, not create arbitrary superposition states. The qubit is the two 
hyperfine levels denoted |0) and |1) in Figured! and to perform operations between 
these states, two-color fields, which are resonant with |0) — > |e) and |1) — > |e) at the 
same time, are employed [9]. Determined by the phase relation, 0, between the two 
fields, two superposition states, shown in Figure Hb, the bright and dark state: 

\B) = ^ 2 (\0)-e^\l)) 

\D) = j= 2 (\0) + e-^\l)) 1 1 

can be defined. Similarly to electromagnetically induced transparency (EIT), the dark 
state does not interact with the light. Thus, the two-color pulse will only target the 
bright state, and two consecutive transfers, first from \B) — > |e) and then from |e) — > \B), 
but going down along a different path as illustrated in Figure |4b, will add a phase angle 
enclosed by the rotation, 9, to the bright state: \B) — > e id \B). This operation can be 
rewritten in the computational basis (|0),|1)), and then becomes 

cos(fl/2) ze^sin(0/2) 1 
ze-^sin(0/2) cos(#/2) ' ^ ' 
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Figure 5: (color online) a) shows the projective measurements on the x-, y- and z-axis 
respectively. The two left peaks are both transitions from the |0) state and correspond 
to a +1 measurement value, while the two right peaks similarly are transitions from the 
|1) state and correspond to a —1 value with the respect to the current axis. Panel b) 
shows the result of a quantum state tomography of the six states along the positive and 
negative axes. 



which is a matrix describing an arbitrary rotation around any vector on the equator of 
the Bloch sphere. By picking suitable phase angles, one can thus perform any qubit 
gate operation [9]. 

These dark state gates were implemented and the qubit superposition states that 
were created as a result of the rotations, were characterized using quantum state 
tomography In Figure [5^ the readout of a particular superposition state, ^= (|0) — |1)) 
is shown as an example. The three curves correspond to projection measurements on 
the x-, y- and z-axis respectively. As expected, we see that both the y and the z 
measurement contain equal contributions from either state, but the x measurement 
gives a projection onto one of the states. States corresponding to the six different 
(positive and negative) axes on the Bloch sphere, were prepared, and characterized, and 
the resulting state vectors inscribed in a Bloch sphere are displayed in Figure [5b. 

The z-projection was done through an absorption measurement, where the two 
leftmost peaks in Figure [5] correspond to transitions from the |0) state and the two 
rightmost peaks correspond to transitions from the |1) state. In order to do the x- and 
^-measurements, suitable rotations were first performed to rotate those axes onto the 
z-axis, such that an absorption measurement now gives the same result as a direct x- or 
y-measurement would have, in the original basis. From those measurements the density 
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matrix can then be calculated through (see e.g. [TP"] ) 

_ tr(p)I + tvjXpjX + ti(Yp)Y + tr(Zp)Z 
P ~ 2 ' [6) 

where X, Y and Z correspond to the Pauli matrices, which together with the identity, 
/, span the 2x2 qubit space in the density matrix representation. The above equation 
is useful because the traces correspond to actual physical measurements. For example, 
tr(Zp) means measuring the projection of the unknown state onto the Z axis, which 
will yield a value between -1 and +1, and the same for the other axes. The fidelity of 
the entire procedure can then be calculated from 

Ftot — i^theor \pexp \ ^theor) > (4) 

i.e. the overlap between the theoretical and the experimental representations. For the 
six states displayed in Figure [5b, the total fidelities were all between 0.84 and 0.92. As 
mentioned above the total procedure consists of two rotations, first one to create the 
state, then another one to perform the measurement in the right basis. Thus, the fidelity 
for a single qubit gate can be expressed as F gate = y/F tot , and we see that the qubit 
gate fidelities are in the range of 0.91 < F gate < 0.96. More information and discussion 
regarding these results can be found in [HI [IH] • 



5.2. Optimal control pulses 

Even though the dark state qubit gates were realized with good fidelity there are still 
reasons to search for further improvements. The total duration of the dark state pulses 
that performs one arbitrary qubit rotation is close to 18 |xs. The employed laser power 
on the other hand, typically provides Rabi frequencies in the order of 1 MHz, which 
indicates that the qubit rotations could probably be performed an order of magnitude 
faster, with an optimal scheme. To reduce the operation time, pulses obtained through 
optimal control theory calculations are currently investigated. 

Optimal control theory provides a framework for finding the best set of controls to 
steer a system so that a desired target state or unitary gate is implemented. In some 
special cases these controls can be determined analytically, see for example [201 12T1 122]. 
In addition, powerful numerical methods [23] are available that make it possible to 
explore the physical limits of time-optimal control experiments in cases where no 
analytical solutions are known. These methods have so far been successfully applied 
to spin systems [2U [251 [2E] and superconducting qubits [27], and trapped ions [28] . 
Optimized pulses can be designed to account for experimental errors, and, if a realistic 
model is available, can include multiple system levels and transitions. 

The basic idea of optimal control algorithms is to minimize a cost function, which is 
determined by the physical system at hand. In the present case, this is given by the light- 
field interaction with the Pr 3+ ion, and in addition to this, a number of experimental 
limitations, such as AOM maximum bandwidth or laser power fluctuations, are enforced 
as well. The algorithm then finds a pulse shape that performs a particular operation to a 
specified minimum fidelity, while staying within the restrictions. It can be noted that the 



Extracting high fidelity quantum computer hardware from random systems 



10 



(a) amplitude envelope (b) phase 




time (ijs) time (/js) 



Figure 6: (color online) Two pulse shapes with different bandwidths, both designed by 
optimal control algorithms to perform efficient state-to-state transfers. The left side 
shows the detected beating and the amplitude envelope, while the right side shows the 
phase evolution. In both cases the red line is the theoretically desired shape, while the 
blue line is the experimentally detected shape, and in both cases they are matching very 
well. Note that in the phase plot, the phase does not have any meaning outside the 
pulse envelope (since the amplitude there is zero). 



solution is not necessarily unique, or even at the global maximum, and the shapes that 
comes out of the calculations are then often strange-looking, without obvious intuitive 
explanations. In order to characterize how well the AOM system could reproduce these 
strange shapes, an interference experiment was set up, where the AOM pulse-shaping 
system is put in one arm and the light passing through this part then interferes with 
the light from an undisturbed part of the laser. The result is a beating, from which 
both the amplitude envelope and the phase evolution of the pulses can be obtained via 
a Fourier transform. Figure [6] shows two examples of such pulses. In panel a) of this 
figure, the green, oscillating curve is the detected beating, the blue line is the amplitude 
envelope of the beating, which is almost perfectly overlapping with the red line, which 
is the theoretically desired shape (may be hard to see if not in color). 

The upper and lower pulse differs in allowed bandwidth. The lower pulse has a 
bandwidth of about 16 MHz, while the upper pulse has been restricted to work with 
a bandwidth of only 2 MHz. Ideally, one would want as high bandwidth as possible, 
since this enables pulses with shorter duration. The right side of Figure [6] displays the 
phase evolution of the detected pulses as well as the desired phase chirp, and again, 
there is almost perfect overlap. In order to test the optimal control results, a series 
of five pulses were designed, each one with the same objective, to do a state-to-state 
transfer from |0) to |e), but with different bandwidths. The resulting efficiency as a 
function of bandwidth is shown in Figure As can be seen, the efficiency is very good 



Extracting high fidelity quantum computer hardware from random systems 



11 



100 



*- 





h- 20 







10° 



10 1 

Bandwidth (MHz) 



Figure 7: Transfer efficiency of |0) — > \e) by means of optimal control pulses, as a 
function of bandwidth. 

for all bandwidths except the highest (16 MHz). The reason for this is not yet fully 
known, but is believed to be related to complications that arise when multiple levels 
are involved. From the level diagram (Figure H]) we see that any pulse that centers 
around the |0) — > |e) transition and has a bandwidth less than 2-4.6 = 9.2 MHz will 
only involve those two levels and nothing else, but any pulse with a higher bandwidth 
will automatically involve at least the two other excited states as well, which creates 
a much more complicated situation. For two-level schemes, other groups have also 
obtained very good results, for example in ion traps [H] . But clearly, there should also 
be suitable optimal pulses for three-level configurations. Thus, the good results for the 
simpler pulses together with the potential of an order of magnitude faster pulses make 
the optimal control technique seem promising for the future. 

6. Single-ion qubits 

6.1. Scalability of the ensemble qubit approach 

Scaling to large number of qubits is a key issue for a quantum computer. In the present 
system, the qubit-qubit operation is achieved by the dipole-dipole interaction between 
two neighboring ions. As stated in Section [2] they must physically sit close enough to 
each other so that when the state of one ion is changed, e.g. from the ground state 
to the excited state and the electric field from the static dipole moment changes, this 
will shift the neighboring ion out of resonance with its original transition frequency. In 
this way, closely lying ions can control each other. For an ion absorbing at a selected 
frequency, Ui, there will be several ions in the vicinity that can control this ion. In the 
ensemble approach, for two qubits, % and j, at transition frequencies z/j and Uj, there is 
a probability, p, that for an arbitrary ion, a;, in qubit i, there will be an ion in qubit j 
that can control Oj. For n qubits, if iV denotes the number of ions in a qubit that can 
be controlled by ions in the other n — 1 qubits, we get that N oc p n ~ x . For the materials 
studied so far p is in the order of 1%, which means that for a 5-qubit quantum register 
there will only be 1 out of 10 s ions that are useful in each qubit. 

Several remedies were pointed out to improve the scaling property by pushing p > 1, 
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e.g. choosing ions with large static dipole moment, increasing the dopant concentration 
[29] . constructing the qubit-qubit interaction by using a mediator bus ion [12] or adding 
a specially selected readout ion to detect the state of single ions, thereby enabling the use 
of single ions instead of an ensemble of ions for each qubit. Considering both advantages 
and disadvantages in each scheme the single ion readout approach is believed to be the 
easiest way to achieve the scalability. 

6.2. The readout-ion approach 

Inside any small volume within a rare-earth-ion doped crystal there is a high probability 
to find several strongly interacting ions which can control each other and each of them 
can represent one qubit. A readout scheme has been suggested for detecting the qubit 
state in such a single ion qubit system through the dipole-dipole interaction [12]. The 
dopant concentration of the readout ion should be very low to make sure there is on 
average only one readout ion interacting with the laser field. To fulfill its role, the 
readout ion is supposed to meet several requirements: (i) The absorption spectrum 
should be well separated from that of the qubit ion so that the readout procedure does 
not affect the qubit state, (ii) The excited state lifetime should be short and there should 
not be any trapping state. Thus, the excitation may be cycled many times producing 
many photons and giving a strong detection signal, (iii) The homogeneous linewidth 
should be narrow so that when a nearby qubit ion is transferred between the ground 
and excited states, the dipole-dipole interaction is able to shift the readout ion in or out 
of resonance with the readout beam frequency, which consequently turns on or off the 
fluorescence, (iv) The inhomogeneous linewidth should be large which provides a larger 
number of frequency channels. 

6.3. Find and search schemes 

Provided the spectral information of the readout ion is characterized, the nearest qubit, 
ql, and further qubit chain structures in the vicinity of the readout ion, can be mapped 
out using the fluorescence signal from the readout ion. After detecting the readout 
ion fluorescence at a certain excitation frequency, u r , one would keep the readout laser 
working at u r , while scanning the qubit laser across the inhomogeneous qubit transition 
frequency to find a frequency, ui, where the fluorescence signal stops. This would 
mean that a qubit ion with the addressing frequency V\ can control the readout ion, 
as illustrated in Figure [SJ With the qubit and readout laser working at v\ and v r 
respectively, the presence or absence of fluorescence tells whether ql is in the |1) or 
the |0) state. The general case is that the qubit is in a superposition state, then one 
has to do the quantum state tomography measurements as described in Section 15.11 
then execute the above readout procedure to characterize the probability distribution 
between the |0) and |1) state for each projection. After the excitation frequency of ql is 
known one can then extend the readout scheme to find out which qubit, q2, can control 
ql through the qubit-qubit interaction and further, consecutively map out the whole 
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Figure 8: (color online) The dipole-dipole interaction between a qubit ion and a near- 
lying readout ion. If the qubit ion ql is excited from the |0) state to the |e) state at 
frequency v\, the electric field generated by the dipole will shift the transition frequency 
of the readout ion by an amount 8u, after which excitation at frequency v r will have no 
effect and the readout ion fluorescence will be turned off. 



interacting qubit chain and read out all qubits. 
6.4- Investigation for a readout ion 

So far Ce 3+ is believed to be the most promising candidate for the readout ion. The 
excitation wavelength for the 4f — 5d zero-phonon-line (ZPL) of Ce 3+ ion randomly 
doped in the YSO-crystal is around 371 nm, which is well separated from the qubit ion 
transition frequency 606 nm (Pr) or 580 nm (Eu). The Ce 3+ homogeneous linewidth 
is determined by detecting the sum frequency signal from a saturation spectroscopy 
experiment with different amplitude modulation frequencies on the pump and probe 
beams. The measured homogeneous linewidth was as narrow as 3 MHz. This is a 
very encouraging result, because based on a reasonable assumption about the Ce dipole 
moment, the Pr ion should be able to shift the Ce transition line by 30 MHz when they 
are separated by 7 nm [30]. With a 50 ns lifetime [3U [32] this also means that the 
homogeneous linewidth is lifetime-limited which is the ideal case for this readout ion 
concept. 

The inhomogeneous linewidth of the ZPL in a crystal with 0.08% cerium 
concentration is measured to be 80GHz at 4K. However, the fluorescence at the center 
of the ZPL is only twice the background signal outside the ZPL. The exact origin of this 
background is unknown, but it could for example be phonon sideband transitions, which 
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are stronger for 4f — 5d transitions than for 4f — 4f transitions. Clearly this will affect the 
readout fidelity. Assume the detector quantum efficiency is 10% and the fluorescence 
collection efficiency is 30% [33], and a 7r-pulse is applied on the |0) — > |e) qubit transition. 
As the readout laser is turned on, there will then be about 100 fluorescence photons 
detected during the optical lifetime (150 \is) of the Pr 3+ ion, if the qubit was in its 
|1) state. On the other hand, if the Pr 3+ ion was in its |0) state, the ZPL background 
fluorescence might still give around 50 photons due to off-resonant excitations (the 
background signal). However, the difference in the number of photons is certainly still 
sufficient for separating the two cases. One way to reduce the background signal, would 
be to decrease the number of Ce ions within the laser focus. 

7. Scalability 

When the single-ion readout system is scaled up by adding new qubits to the (possibly 
branched) qubit ion chain there will, although there are many frequency channels 
available within the inhomogeneously broadened qubit transition line, eventually be 
a new qubit ion (Q new ) that have the same frequency as one of the qubit ions (Q id) 
already in the chain. In this situation, the two ions (Qdd and Q nC w) cannot any longer 
be individually addressed in frequency space. It is conceivable that there are algorithmic 
solutions to this problem, where single qubit operations could be replaced by two-qubit 
operations, so that operations are applied to "the qubit with frequency ^, sitting next 
to a qubit with frequency vf . As an alternative, the problem of coinciding transition 
frequencies can be removed by mounting (closely spaced) electrodes onto the crystal. 
The part of the qubit chain on which the operations are carried out, is selected by 
applying a voltage on the appropriate electrode. Eu in Y 2 SiOs has a Stark coefficient of 
35 kHz/(V/cm) [31], and ions close to the electrode can be shifted into resonance with 
the excitation pulse. A 1 MV/cm field would shift the transition frequency 35 GHz, a 
detuning much larger than the inhomogeneous transition line width. Simulations show 
that 20 nm long electrodes separated by 40 nm on a surface may give shifts of the order 
of 30 GHz within a 60 nm region while ions 20 nm further away would shift less than 10 
GHz [12]. However, still smaller electrodes and electrode separations would be favorable 
and while this is in line with the current technological development we will also now 
analyze the possibility of entangling spatially remote few-qubit systems. 

Thus, we consider spatially separated qubit chains, each containing a not too large 
number of qubits. There are several schemes for entangling spatially remote few-qubit 
systems. To optically entangle a qubit in one system with a qubit in another system 
the key point is to configure the photon detection (or absorption) process such that 
when a photon is detected (absorbed) it cannot in principle be determined which of the 
two systems that emitted (absorbed) the photon [351 [36]. If t ne set-up is appropriately 
configured, the photon detection/absorption entangles the two qubits. 

Scalable quantum computing is then achievable also with errors in the remote 
entanglement in the order of several percent, provided that the error probability for 
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operation within each of the qubit chain is below 10~ 4 [H7J. Remote entanglement 
processes often correspond to a probabilistic event, which is repeated until it succeeds. 
Clearly, schemes with very high probability to succeed in a single try are preferable. 
Such schemes can for example be devised using optical micro- cavities [38], e.g., toroidal 
micro-resonators [39j HQ] and can also be improved by the use of cluster states [H]. 

We consider a rare earth crystal with a Ce readout ion, close to the surface facing 
a toroidal microcavity. Around this readout ion there will be a number of qubit ions 
that interact and can be read out, as described in Section EJ One of these qubit ions is 
selected as the entangler ion. The purpose is now to entangle entangler ions in different 
rare-earth crystals. Such a scheme requires that also rare-earth ions close to the crystal 
surface retain their favorable coherence properties, and preliminary data indicate that 
this indeed is the case |42j. 

Consider Figure El the interaction between the cavity and the ion is adjusted such 
that a single photon corresponds to a 7r-pulse [43]. Using tapered fibres, a single photon 
can be coupled into the toroidal microcavity with an efficiency (ideality) > 99.98% 
[H]. The entangler ion in the first crystal, A, is prepared in the superposition state 
(|0) + \1))a and the entangler ion in the second crystal, B, is prepared in state |1)b- 
As experimentally demonstrated [15] , the ion transition frequencies to an excited state, 
|e), can readily be tuned by an electric field. Thus, by tuning the |0) — > |e) transition 
for the ion in A to the same frequency as the |1) — > |e) transition for the ion in B using 
electric fields (see right hand side of Figure [9]) and sending a photon in the tapered fibre 
at this frequency in the direction indicated in the figure, the two ions will be prepared 
in the state i\e) A |1) B + |1) A \e) B ). Applying local 7r-pulses on each of \e) A —> \0) A and 
\e) B — > |0) at the respective input ports for the qubit and the readout ion, then creates 
the Bell state |01) + |10). 

In a more complete version of the approach in Figure [91 there would be a large 
number of microresonators (with their corresponding crystals) along the tapered fibre. 
Which two crystals that are selected for any specific entanglement process can be 
controlled by e.g. adjusting the distance between the tapered fibre and the resonator 
(cavity) using e.g. micro-actuators. Alternatively, the cavity can be tuned off resonance 
with the transition and photon frequency. Tuning the cavity can be carried out using 
expansion due to temperature. For example, electrodes can be connected to the toroid 
and the substrate. The temperature and thereby also the frequency change will then be 
proportional to the dissipated power [16] . In general the system could be quite flexible 
and versatile, as an example, ions can be tuned in and out of resonance with the cavity 
using an electric field across the crystal. 

As earlier stated, the qubit coherence times can be 30 seconds for Pr ions [7J and 
for Eu ions coherence times of the order of an hour have been predicted [17] . A strength 
with the scheme outlined above is therefore that there is ample time both to create 
entanglement and to carry out operations once the systems have been entangled. 

To compare with some existing data, free Cs atoms can have a single photon 
Rabi frequency of 100 MHz for Cs atoms 45 nm from the toroidal surface [H]. For 
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Figure 9: (color online) Schematic view of an approach to entangle spatially separate 
few qubit processors. 

p r 3+ : Y 2 gj0 5 which has an oscillator strength almost six orders of magnitude lower than 
free Cs atoms [M], would then yield a Rabi frequency of about 150 kHz. This value is 
too low for a single photon to act as a 7r-pulse and we are therefore looking at techniques 
for enhancing the interaction between an ion just below the surface and an evanescent 
field. 

8. Concluding remarks 

To summarize, single-qubit operation fidelities above 90% are obtained for ensemble 
qubits and > 99% single qubit operation fidelity should be readily obtained with single- 
ion qubits. However, the single ion qubit scheme is contingent upon the ability to read 
out the hyperfine state of single rare-earth ions and work to develop this capacity is still 
ongoing. Although, based on current data in the investigation of using Ce 3+ as a readout 
ion for such a scheme, there is reason to be quite optimistic about the possibilities to 
succeed. Finally, it may be pointed out that while quantum state storage in rare-earth- 
ion-doped materials is pursued vigorously by several groups and now demonstrates 
excellent progress [19] including storage efficiencies above 30% [50], only two groups 
[5TT E] have published experimental efforts on quantum computing in rare earth crystals 
and considering the limited efforts, substantial progress has been made. 



Extracting high fidelity quantum computer hardware from random systems 



17 



Acknowledgements 

This work was supported by the Swedish Research Council, the Knut and Alice Wallen- 
berg Foundation, the European Commission through the integrated project QAP and 
the Crafoord Foundation. 

[1] Y. Nakamura, Y. A. Pashkin, and J. S. Tsai. Coherent control of macroscopic quantum states in 

a single- cooper-pair box. Nature, 398(6730):786-788, 1999. 
[2] J. H. Plantenberg, P. C. de Groot, C. J. P. M. Harmans, and J. E. Mooij. Demonstration of 

controllcd-not quantum gates on a pair of superconducting quantum bits. Nature, 447:836, 

2007. 

[3] D. Loss and D. P. DiVinccnzo. Quantum computation with quantum dots. Phys. Rev. A, 
57(1):120-126, 1998. 

[4] B. E. Kane. A silicon-based nuclear spin quantum computer. Nature, 393(6681):133-137, 1998. 
[5] F. Jelezko, T. Gaebel, I. Popa, M. Domhan, A. Gruber, and J. Wrachtrup. Observation of coherent 

oscillation of a single nuclear spin and realization of a two-qubit conditional quantum gate. Phys. 

Rev. L., 93:130501, 2004. 

[6] N. Ohlsson, R. K. Mohan, and S. Kroll. Quantum computer hardware based on rare-earth-ion- 
doped inorganic crystals. Opt. Commun., 201(l-3):71-77, 2002. 
[7] E. Fraval, M. J. Scllars, and J. J. Longdcll. Dynamic decoherence control of a solid-state nuclear- 

quadrupolc qubit. Phys. Rev. Lett., 95(3):030506, July 2005. 
[8] L. Rippc, B. Julsgaard, A. Walther, Yan Ying, and S. Kroll. Experimental quantum state 

tomography of a solid state qubit. Phys. Rev. A, 77:022307, 2008. 
[9] I. Roos and K. Maimer. Quantum computing with an inhomogcncously broadened ensemble of 

ions: Suppression of errors from detuning variations by specially adapted pulses and coherent 

population trapping. Phys. Rev. A, 69(2):022321, 2004. 
[10] L. Rippe, M. Nilsson, R. Klicber, D. Suter, and S. Kroll. Experimental demonstration of efficient 

and selective population transfer and qubit distillation in a rare-earth-metal-ion-doped crystal. 

Phys. Rev. A, 71(6):062328, June 2005. 
[11] N. Timoncy, V. Elman, S. Glascr, C. Weiss, M. Johanning, W. Neuhauser, and C. Wundcrlich. 

Error-resistant single-qubit gates with trapped ions. Phys. Rev. A, 77:052334, 2008. 
[12] J. H. Wesenberg, K. Molmer, L. Rippc, and S. Kroll. Scalable designs for quantum computing 

with rare-earth-ion-doped crystals. Physical Review A, 75(1):012304, January 2007. 
[13] R. W. Equall, R. L. Cone, and R. M. Macfarlane. Homogeneous broadening and hyperfine structure 

of optical transitions in Pr 3+ : Y 2 Si0 5 . Phys. Rev. B, 52:3963-3969, 1995. 
[14] Y. Sun, C. W. Thicl, R. L. Cone, R. W. Equall, and R. L. Hutchcson. Recent progress in developing 

new rare earth materials for hole burning and coherent transient applications. J.Lumin., 98(1- 

4):281-287, 2002. 

[15] E. Fraval, M. J. Scllars, and J. J. Longdcll. Method of extending hyperfine coherence times in 
Pr 3+ : Y 2 Si0 5 . Phys. Rev. Lett, 92(7):077601, 2004. 

[16] B Julsgaard, L Rippe, A Walther, and S Kroll. Understanding laser stabilization using spectral 
hole burning. Optics Express, 15:11444-11465, 2007. 

[17] O. Guillot-Noel, Ph. Goldncr, F. Beaudoux, Y. Le Du, J. Lejay, A. Amari, A. Walther, L. Rippc, 
and S. Kroll. Hyperfine structure and hyperfine coherent properties of praseodymium in single- 
crystalline La2{WOi)3 by hole-burning and photon-echo techniques. Phys. Rev. B, 79:155119, 
2009. 

[18] Andreas Walther. Coherent processes in rare- earth-ion- doped solids. PhD thesis, Lund University, 
2009. http://www-atom.fysik.lth.se/QI/publications/Theses/Thesis-AW.pdf. 

[19] M. A. Nielsen and I. L. Chuang. Quantum computation and quantum information. Cambridge 
University Press, United Kingdom, 2000. Eq. 8.148, Chapter 8.4.2. 



Extracting high fidelity quantum computer hardware from random systems 



18 



[20] Navin Khaneja, Roger Brockett, and Steffen J. Glaser. Time optimal control in spin systems. 

Phys. Rev. A, 63(3):032308, Feb 2001. 
[21] Navin Khaneja, Burkhard Luy, and Steffen J. Glaser. Boundary of quantum evolution under 

dccoherence. Proc. Natl. Acad. Sci. USA, 100:13162-13166, 2003. 
[22] Robert Fisher, Haidong Yuan, Andreas Sporl, and Steffen Glaser. Time-optimal generation of 

cluster states. Physical Review A (Atomic, Molecular, and Optical Physics), 79(4):042304, 2009. 
[23] Navin Khaneja, Timo Reiss, Cindic Kchlet, Thomas Schulte-Hcrbrggcn, and Steffen J. Glaser. 

Optimal control of coupled spin dynamics: design of nmr pulse sequences by gradient ascent 

algorithms. Journal of Magnetic Resonance, 172(2):296 - 305, 2005. 
[24] Thomas E. Skinner, Timo O. Reiss, Burkhard Luy, Navin Khaneja, and Steffen J. Glaser. 

Application of optimal control theory to the design of broadband excitation pulses for high- 
resolution nmr. Journal of Magnetic Resonance, 163(1):8 - 15, 2003. 
[25] Kyryl Kobzar, Thomas E. Skinner, Navin Khaneja, Steffen J. Glaser, and Burkhard Luy. Exploring 

the limits of broadband excitation and inversion pulses. Journal of Magnetic Resonance, 

170(2):236 - 243, 2004. 

[26] Kyryl Kobzar, Burkhard Luy, Navin Khaneja, and Steffen J. Glaser. Pattern pulses: design of 

arbitrary excitation profiles as a function of pulse amplitude and offset. Journal of Magnetic 

Resonance, 173(2):229 - 235, 2005. 
[27] A. Sporl, T. Schulte-Herbriiggen, S. J. Glaser, V. Bergholm, M. J. Storcz, J. Ferber, and F. K. 

Wilhelm. Optimal control of coupled josephson qubits. Physical Review A (Atomic, Molecular, 

and Optical Physics), 75(1):012302, 2007. 
[28] V. Nebendahl, H. Haffner, and C. F. Roos. Optimal control of entangling operations for 

trapped-ion quantum computing. Physical Review A (Atomic, Molecular, and Optical Physics), 

79(1):012312, 2009. 

[29] Mattias Nilsson. Coherent interactions in rare-earth-ion-doped crystals for applications in quantum 
information science. PhD thesis, Division of Atomic Physics, LTH, 2004. 

[30] Christian Ahlholm. Simulation of the single ion readout concept. Master's thesis, Lund University, 
2008. Lund Reports on Atomic Physics, LRAP-402. 

[31] Julio E. Hernandez. Spectroscopy in Ce 3+ : Y2S1O5. a preliminary investigation for a single ion 
readout scheme for quantum computation with rare-earth ion doped crystals. Master's thesis, 
Lund University, 2006. Lund Reports on Atomic Physics, LRAP-356. 

[32] Tuomas Aitasalo, Jorma His, Mika Lastusaari, Janina Legendziewicz, Janne Niittykoski, and 
Fabienne Pell. Delayed luminescence of Ce 3+ doped Y^SiO^. Optical Materials, 26(2):107 - 
112, 2004. Proceedings of the Second International Conference on Sol-Gel Materials. Research, 
Technology, Applications (SGM'03). 

[33] Viktor Palm. Compact simplified design of optical cell for single-molecule spectroscopy. Rev. Sci. 
Instrum., 70:2957, 1999. 

[34] F. R. Graf, A. Renn, G. Zumofen, and U. P. Wild. Photon-echo attenuation by dynamical processes 
in rare-earth- ion-doped crystals. Phys. Rev. B, 58(9):5462-5478, 1998. 

[35] L.-M. Duan and H. J. Kimble. Efficienct engineering of multiatom entanglement through single- 
photon detections. Phys. Rev. Lett, 90:253601, 2003. 

[36] S. D. Barett and P. Kok. Efficient high-fidelity quantum computation using matter qubits and 
linear optics. Phys. Rev. A, 71:060310, 2005. 

[37] L. Jiang, J. M. Taylor, A. S. S0renssen, and M. D. Lukin. Distributed quantum computation 
based on small quantum registers. Phys. Rev. A, 76:062323, 2007. 

[38] J. I. Cirac, P. Zoller, H. J. Kimble, and H. Mabuchi. Quantum state transfer and entanglement 
distribution among distant nodes in a quantum network. Phys. Rev. Lett., 78(16):3221-3224, 
1997. 

[39] K. J. Vahala. Optical microcavitics. Nature, 424:839, 2003. 

[40] S. M. Spillanc, T. J. Kippcnberg, K. J. Vahala, K. W. Goh, E. Wilcut, and H. J. Kimble. Ultrahigh- 
q toroidal microrcsonators for cavity quantum electrodynamics. Phys. Rev. A, 71:013817, 2005. 



Extracting high fidelity quantum computer hardware from random systems 19 

[41] R. Raussendorf and H. J. Briegel. A one-way quantum computer. Phys. Rev. L., 86:5188, 2001. 

[42] J. Longdell et al. Cavity QED with rare earths. 2008. Workshop in Bozcman, Montana. 

[43] Y.-F. Xiao, Z.-F. Han, Y. Yang, and G.-C. Gou. Quantum cpf gates between rare earth ions 

through measurement. Phys. Lett. A, 330:137, 2004. 
[44] S. M. Spillane, T. J. Kippcnberg, O. J. Painter, and K. J. Vahala. Ideality in a fiber-taper- 

couplcd microresonator system for application to cavity quantum electrodynamics. Phys. Rev. 

Lett, 91(4):043902, Jul 2003. 
[45] R. M. Macfarlane and A. J. Meixner. Electric-ficld-modulatcd photon echoes in Pr 3+ : YAIO3. 

Opt. Lett, 19:987, 1994. 

[46] D. Armani, B. Min, A. Martin, and K. J. Vahala. Electrical thermo-optic tuning of ultrahigh-q 

microtoroid resonators. Appl. Phys. Lett., 85:5439, 2004. 
[47] J. J. Longdell, A. L. Alexander, and M. J. Scllars. Characterization of the hyperfme interaction 

in europium-doped yttrium orthosilicate and europium chloride hexahydrate. Phys. Rev. B, 

74:195101, 2006. 

[48] T. Aoki, B. Dayan, E. Wilcut, W. P. Bowen, A. S. Parkins, T. J. Kippenberg, K. J. Vahala, and 
Kimble H. J. Observation of strong coupling between one atom and a monolithic microresonator. 
Nature, 443:671, 2006. 

[49] H. Riedmatten, M. Afzelius, M. U. Staudt, C. Simon, and N. Gisin. A solid-state light-matter 

interface at the single-photon level. Nature, 456:773, 2008. 
[50] A. Amari, A. Walther, M. Sabooni, M. Huang, S. Kroll, M. Afzelius, I. Usmani, H. de Riedmatten, 

and N. Gisin. Towards an efficient atomic frequency comb quantum memory. arXiv:0911.2145 

[quant-ph], 2009. 

[51] J. J. Longdell and M. J. Sellars. Experimental demonstration of quantum-state tomography 
and qubit-qubit interactions for rare-earth-metal-ion-based solid-state qubits. Phys. Rev. A, 
69(3):032307, 2004. 



